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Figure  1.  Sample  strip  chart  of  raw  data  output  in  the  field  for  equipment  diagnosis  and  4 
quick-look  data  assessment.  Panels,  from  top  to  bottom,  show  differential  carrier 
phase  (DCP)  between  LI  (1575.4  MHz)  and  L2  (1227.6  MHz),  differential  group 
delay  (DGD)  between  the  same  two  signals,  and  their  received  intensities  (P2  and 
Pj  for  L2  and  LI,  respectively)  in  dB  relative  approximately  to  their  means  over 
the  26  minutes  shown.  Note  the  inversion  (0  on  top)  of  the  DGD  scale.  Data 
recorded  at  Lerwick  on  31  December  1989. 

Figure  2.  Example  of  spreadsheet  catalogue  and  quick-look  synopsis  of  TEC  and  scintil-  5 
lation  activity,  showing  data  for  the  month  in  which  the  record  shown  in  Fig.  1 
was  collected. 

Figure  3.  Intensity  scintillation  indices,  S4,  computed  for  L2  (top)  and  LI  (bottom)  8 

from  the  raw  data  illustrated,  in  part,  in  the  bottom  two  panels  of  Fig.  1.  Time 
span  is  13:13  through  1 6: 13  UT.  The  three  spikes  in  S4  correspond  to  the  three 
patches  of  scintillation  in  Fig.  1.  The  large  spike  in  the  L2  value  occurred  later 
than  the  end  of  Fig.  I  and  is  an  obvious  artifact. 

Figure  4.  TEC  values  corresponding  to  S4  records  shown  in  Fig.  3,  derived  from  the  9 

DGD  (smoother  of  the  two  overlapping  traces  at  the  top)  and  DCP  (more  finely 
structured  trace)  measurements  shown,  in  part,  in  Fig.  1.  Lower  solid  trace  is 
equivalent  vertical  TEC  derived  from  DCP  by  removing  effect  of  changing  slant 
path-length,  x  and  +  indicate  satellite  azimuth  and  elevation,  respectively. 

Time  span  is  13:14  through  16:14  UT. 

Figure  5.  Synoptic  chart  of  equivalent  vertical  TEC  values  (solid  curves)  derived  from  10 
DCP  and  calibrated  by  means  of  DGD  from  measurements  at  Lerwick  over  the 
24-hour  period  (31  December  1989,  UT)  in  which  the  satellite  track  producing 
the  data  in  Figs.  1 , 3,  and  4  occurred,  x  and  +  indicate  satellite  latitude  and 
longitude,  respectively. 

Figure  6.  Synoptic  chart  01  equivalent  vertical  'FEC  and  satellite  latitude  and  longitude  1 1 
for  24-hour  period  of  20  October  1989,  UT. 
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Figure  7.  Quick-look  spreadsheet  output  from  Lerwick  fo:  September  1989,  indicating  12 
scintillation  and/or  elevated  TEC  activity  on  UT  dates  16,  18, 19,  22,  and  26,  all 
of  which  were  geomagnetically  disturbed  days. 

Figure  8.  Synoptic  chart  of  equivalent  vertical  TEC  and  satellite  latitude  and  longitude  1 3 
for  24-hour  period  of  22  September  1989,  UT.  Scintillation  occurred  during  UT 
hours  14  and  15,  in  the  presence  of  level-1  TEC  activity. 

Figure  9.  Latitudinal  scans  of  TEC  observed  at  Lerwick  by  means  of  polar-orbiting  14 

"Transit"  satellites  between  09:06  and  09:18  UT  (shallow  curve)  and  between 
15:18  and  15:29  UT  (steeper  curve)  on  22  September  1989.  The  former  was 
approximately  along  the  353-deg  meridian.  (Data  provided  by  L.  Kersley  via 
G.  Bishop.) 


Figure  10.  Synoptic  chart  of  equivalent  vertical  TEC  and  satellite  latitude  and  longitude  15 
for  24-hour  period  of  26  September  1989,  UT. 

Figure  11.  TEC  ineasured  during  the  first  afternoon  GPS  track  in  Fig.  8.  16 

Figure  12.  Intensity  scintillation  indices,  S4,  measured  during  the  first  afternoon  GPS  17 

track  in  Fig.  8. 

Figure  13.  TEC  measured  during  the  second  afternoon  GPS  track  in  Fig.  10.  Note  18 

strong  quasi-pericdic  signature  centered  on  about  17:27  UT. 

Table.  TEC  Variability  Levels.  6 
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PRELIMINARY  INVESTIGATION  OF  LONGITUDINAL  DIFFERENCES  IN  TEC 
AND  SCINTILLATION  AT  TRANSITION  LATITUDES 


I.  INTRODUCTION  AND  OBJECTIVES 

The  F  layer  of  the  earth's  ionosphere  behaves  rather  differently  in  three  latitudinal 
regions.  The  simplest  -  produced  almost  entirely  by  solar  ultraviolet  (uv)  radiation  and  con¬ 
trolled  largely  by  diffusion  -  is  the  mid-latitude  region,  between  about  20  and  50  deg  geo¬ 
magnetic  latitude.  The  equatorial  F  layer,  at  lower  geomagnetic  latitudes,  also  is  produced  pri¬ 
marily  by  solar  uv,  but  it  is  substantially  more  dynamic  and  structured  on  the  night  side  of  the 
earth  than  is  its  mid-latitude  neighbor.  Specifically,  its  plasma  there  is  subject  to  convective 
interchange  instabilities  analogous  to  the  Rayleigh-Taylor  instability,  well-known  in  neutral  fluid 
dynamics.  Such  instabilities  result  in  substantial  spatial  variation  in  j^lasma  density  and  its  path 
integral.  Total  Electron  Content  (TEC).  At  sufficiently  small  scales,  these  structures  produce 
random  fluctuations  -  "scintillations"  -  in  the  phase,  polarization,  and  intensity  of  radio  signals 
received  via  transionospheric  channels. 

Poleward  of  the  mid-latitude  ionospheric  lies  another  structured  and  still  more  dynamic 
region.  The  high-latitude  ionosphere  is  produced,  not  only  by  solar  uv,  but  also  by  energetic 
particles  (primarily  electrons  with  many  kev  of  energy)  that  precipitate  into  it  from  the  magneto¬ 
sphere.  In  addition  to  providing  a  "storage  bin"  for  these  ptirticles  (also  ultimately  of  solar 
origin),  the  magnetosphere  --  and  its  interaction  with  the  streaming  plasma  and  imbedded  mag 
nctic  field  of  the  solar  wind  -  largely  control  the  high-latitude  F  layer  dynamically.  Not  only  do 
plasma-interchange  instabilities  produce  density  and  TEC  structures  there,  but  the  prirnmy  force 
that  drives  them  -  the  E  x  B  force,  where  E  and  B  respectively  denote  the  local  electric  and 
magnetic  fields  -  dictates  convective  motion  of  the  F-layer  plasma  and  its  structures. 

The  high-latitude  ionosphere  may  be  divided  further  into  the  auroral  and  polar  sub- 
regions,  dominated  respectively  by  particle  and  convective  effects.  Our  interest  in  this  report  is 
in  plasma  stnictures  and  dynamics  in  the  transition  region  between  the  mid-latitude  and  auroral 
ionospheres.  Generally,  the  transition  region  lies  between  about  50  and  65  deg  geomagnetic 
latitude,  but  the  definition  is  not  a  static  one.  By  transition  latitudes,  we  mean  the  region  of  the 
ionosphere  extending  from  the  plasmapause  (the  poleward  tenninus  of  the  mid-latitude  iono¬ 
sphere)  through  the  scintillation  boundary  that  lies  near  the  equatorward  edge  of  the  auroral  oval 
(the  region  of  diffuse  aurora).  It  includes  the  main  plasma  trough  that  separates  the  uv-derived, 
diffusively  driven  mid-latitude  F  layer  and  the  particle-influenced,  convecting  high-latitude  F 
layer.  It  lies  more  equatorward  on  the  night  side  of  the  earth  than  on  the  day  side  and  migrates 
progressively  equatorward  with  increasing  geomagnetic  disturbance.  In  the  North  Atlantic,  the 
transition  region  lies  between  geographic  latitudes  of  about  45  and  70  deg,  somewhat  more 
southerly  on  the  American  side  and  northerly  in  the  European  sector. 
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Northwest  Research  Associates  (NWRA)  has  conducted  a  preliminary  investigation  of  L- 
band  group  delay,  dispersion,  and  scintillation  produced  by  the  ionosphere  on  opposite  sides  of 
the  North  Atlantic.  The  project  entailed  the  recording  of  differential  group  delay  (DGDj, 
differential  carrier  pha.se  (DCP),  and  amplitude  of  the  two  L-band  signals  received  from  satellites 
of  the  Global  Positioning  System  (GPS),  using  receivers  owned  by  the  Geophysics  Laboratory 
(GL)  of  Air  Force  Systems  Command  (now  an  element  of  Phillips  Laboratory).  Measurements 
were  conducted  at  GL,  located  at  Hanscom  AFB,  MA,  and  at  Lerwick,  in  the  Shetland  Islands, 
UK.  While  about  70  deg  apart  in  longitude  and  almost  20  deg  in  geographic  latitude,  these  two 
stations  lie  closer  to  one  another  in  geomagnetic  latitude,  being  respectively  near  55  and  60  deg 
invariant.  The  DGD  and  DCP  data  were  processed  to  identify  dominant  changes  in  group  delay 
and  dispersion,  respectively.  In  addition,  the  signal-strength  records  were  processed  to  evaluate 
L-band  scintillation  in  temis  of  the  fractional  standard  deviation,  S4,  of  received  signal  intensity 
(power;. 

The  objectives  of  the  p.oject  were  (1)  to  detennine  differences  between  the  parameters 
measured  in  the  two  different  longitude  sectors  and  (2)  to  interpret  observed  differences  in  terms 
of  ionospheric  total  electron  content  (TEC)  and  intcnncdiate-scale  density  structures  near  the 
two  stations.  The  ultimate  aim  of  the  research  is  to  understand  the  TEC  and  structural 
differences  within  the  context  of  high-latitude  ionospheric  and  magnetospheric  dynamics.  From 
the  outset,  however,  the  stud>  was  viewed  as  preliminary.  It  has  been  limited  to  identification  of 
whatever  differences  might  arise  between  the  databases  from  the  two  stations.  Systematically 
categorizing  those  differences  will  require  separation  of  individual  records  in  terms  of  the  time 
and  location  of  their  line-of-sight  intersections  with  the  ionosphere. 

II.  RECEIVER  INSTALLATION  AND  DATA  COLLECTION 

The  first  major  task  performed  under  the  contract  was  to  identify  and  arrange  for  use  of 
an  appropriate  field  site  in  the  UK  and  to  deploy  a  GPS  receiver  there.  Lerwdek  was  chosen 
because  of  its  location  in  the  transition  region  between  gcomagnetically  middle  and  high 
latitudes  (beneath  the  main  ion  trough  and  just  equatorward  of  the  usual  location  of  the  auroral 
oval)  and  becau.se  it  would  facilitate  compari.son  with  ionospheric  data  collected  in  that  region  by 
others  (L.  Kersel>,  private  communication).  Pemiission  to  install  and  operate  the  equipment  at 
the  Lerwick  Observatory  was  obtained  from  the  UK  Meteorological  Office  in  August  1989.  The 
primary  receiver  (a  Stanford  Telecommunications  Model  STEL-5010)  deployed  to  the  Shetlands 
station  previously  had  operated  at  Thule,  Greenland.  It  was  calibrated  and  its  operation 
monitored  for  two  weeks  early  in  the  first  quarter  of  the  project,  and  it  was  shipped  to  the  UK  on 
4  September  1989.  It  arrived  at  Lerwick  on  12  September  and,  after  recalibration,  was  placed  in 
service  on  15  September.  It  his  operated  there  since,  with  some  interruptions  due  to  various 
equipnrnt  and  operational  problems.  A  similar  receiver  (STEL-501())  began  operation  for  this 
project  at  Han.scom  on  25  September  1989,  and  it  continues  in  similar  operation  for  related 
purpo.ses  In  late  March  and  early  April  1990,  NWRA  configured  and  installed  a  backup 
receiver  (a  Texas  Instruments  Model  T1-4I00)  at  Lerwick. 
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Each  STEL  receiver  records  the  intensity  of  the  two  signals  (LI  at  1575.4  MHz,  and  L2 
at  1227.6  MHz)  transmitted  from  a  single  GPS  satellite  during  any  given  observing  period,  along 
with  the  DCP  and  DGD  between  them.  A  satellite  is  doppler  tracked  typically  for  a  few  hours, 
and  then  the  receiver  selects  another  satellite  based  on  the  elevation  of  its  forthcoming  track 
across  the  sky.  The  raw  measurements  are  recorded  at  20  samples  per  sec  (sps),  the  DCP  and 
DGD  data  being  smoothed  and  decimated  later  to  10  samples  per  min  (spm)  for  processing.  The 
TI 4100  receiver  at  Lerwick  is  capable  of  recording  data  from  four  satellites  simultaneously. 

An  example  of  raw  data  from  the  STEL  receiver  at  Lerwick  is  contained  in  Fig.  1,  in  the 
fomi  of  a  strip  chart  run  in  the  field  for  diagnostic  and  quick-look  purposes.  The  two  lower 
panels  show  the  LI  (bottom)  and  L2  signal  strengths,  including  some  intensity  scintillation  -  in 
particular,  three  patches  measuring  +  several  dB  and  lasting  a  few  minutes.  The  top  panel  shows 
the  DCP,  subject  to  a  many-7t  ambiguity.  When  the  2n  crossovers  are  connected,  the  DCP 
provides  a  high-resolution  record  of  relative  TEC,  subject  to  .some  fine-scale  distortion  at  times 
of  scintillation.  The  absolute  value  of  TEC  is  provided,  with  far  coarser  resolution,  by  the  DGD, 
shown  (un  an  inverted  scale)  in  the  remaining  panel  in  TEC  units  of  10*^  el/m^. 


Early  in  the  project,  we  laid  out  and  implemented  a  spreadsheet  routine  to  provide  a 
catalogue  and  quick-look  .synopsis  of  activity  levels  from  Lerwick  and  Hanscom.  An  example  is 
shown  in  Fig.  2,  which  contains  the  spreadsheet  output  from  Lerwick  for  the  month  in  which  the 
record  shown  in  Fig.  1  was  collected  (on  its  last  day).  In  the  main  field  of  Fig.  2  and  of  the 
spreadsheet  outputs  from  other  months,  the  single-digit  entries,  1  through  6,  refer  to  six  levels  of 
TEC  variability,  as  determined  from  the  DCP  records.  (A  0  entry  denotes  an  hour  in  which  no 
data  were  collected.)  This  TEC  activity  scale  originally  was  developed  from  visual  inspection  of 
similar  records  from  Thule,  Greenland.  It  is  quantified  in  the  table  on  page  5,  in  terms  of  TEC 
change  measured  in  TEC  units  per  unit  time.  (See  the  reference*  for  sample  records  that  defined 
the  scale  originally.)  Single  asieri.sks  in  the  spreadsheet  outputs  denote  times  that  intensity 
scintillation  was  detected  on  one  (usually,  the  lower)  oi  both  of  the  two  GPS  frequencies.  The 
threshold  of  detectability  for  scintillation  is  approximately  S4  =  0.05. 
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Lenvick  on  31  December  1989. 


ACTIVITY  REPORT 


PLACE:  SHETLAND  MONTH:  DECEMBER  YEAR:  1989 

FOR 
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COtlHENT'S: 

A:  TAPE  STOP. 

B:  GPS  SYSTEM  OFF  THE  AIR. 

C:  NEW  SATELITE  WINDOWS  INSERTED. 

D:  SOME  DATA  OH  CHART  HOT  RECORDED  OH  TAPE. 
E:  SOME  DATA  BROKEN  UP. 

F:  TAPE  CHANGE. 

TOTAL  HOURS  OF  ACTIVITY  LEVEL  1:300 
TOTAL  HOURS  OF  ACTIVITY  LEVEL  2:59 
TOTAL  HOURS  OF  ACTIVITY  LEVEL  3:16 
TOTAL  HOURS  OF  ACTIVITY  LEVEL  A: 12 
TOTAL  HOURS  OF  ACTIVITY  LEVEL  5:0 
TOTAL  HOURS  OF  ACIIVITY  LEVEL  6:0 


ADDITIONAL  COMMENT'S. 
•  =  L1/L2  SCINTILLATION. 

•*=  NO  ABSOLUTE  TEC. 

T  =  TAPEDRIVE  PROBLEM. 

H  =  HPIB  IHIERFASE  PROBLEM. 

V  =  FLOPPY  OISK/DRIVE  PROBLEM. 

W  =  WRONG  SATELITE  WINDOWS. 

P  =  POWER  FAILURE. 

Q  =  ANTENNA  PROBLEM. 


TOTAL  HOURS  OF  ACTIVITY  LEVEL  2  AND  ABOVE:  87 
TOTAL  HOURS  OF  ACTIVITY  LEVEL  3  AND  ABOVE:  28 


Figure  2.  Example  of  .spread.slieet  catalogue  and  quick-look  synop.si.s  of  TEC  and  scintillation 
activity,  showing  data  for  the  month  in  which  the  record  shown  in  Fig.  1  was  collected. 
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Table.  TEC  variability  levels. 


Level 

Variabilitv 

Comment 

Quiet 

<  0.5  TEC  u/10  min 

No  detectable  activity. 

Small 

<  0.5  TEC  u/  5  min 

Detectable  TEC  variation, 
detectable  scintillation. 

Moderate 

<  1.0  TEC  u/  5  min 

Smooth  rolling  structures. 

Disturbed 

<  1.0  TEC  u/30  sec 

Sharp,  short  gradients. 

Very 

Disturbed 

>  1.0  TEC  u/30  sec 

Very  sharp,  sustained  gradients. 

Extremely 

Disturbed 

>>  1.0  TEC  u/30  sec 

Extremely  large  gradients. 

The  catalog  resulting  from  the  .spreadshcct.s  shows  that,  from  September  1989  through 
November  1990,  some  3700  hours  of  data  were  collected  at  Lerwick  and  approximately  6790 
hours  of  data  were  collected  at  Hanscom.  The  main  interruption  in  data  continuity  occurred  at 
Lerwick  between  mid-June  and  mid-September  1990.  There  was  no  extended  interruption  at 
Hanscom.  Data  interruptions  at  Lerwick  stemmed  from  several  instrument  problems,  as  well  as 
en  route  delays  in  shipping  of  blank  tapes  to  the  site.  The  former  included  temporary  loss  of  an 
antenna  during  a  severe  windstomi  in  the  Shetland  Islands.  Some  problems  did  not  result  in  total 
loss  of  data,  since  some  data  were  recorded  on  strip  charts  and  floppy  disks,  as  well  as  on  the 
primary  magnetic-tape  medium. 

Certain  problems  degraded  data  in  such  a  way  that  the  degradation  could  be  overcome  in 
subsequent  processing.  One  example  involved  calibration  of  the  DCP  records  by  means  of  DGD 
to  permit  measurement  of  absolute  TEC.  The  DGD  measurement  is  obtained  from  the  relative 
time  delay  between  coded  messages  (the  so-called  P  code)  telemetered  on  the  two  frequencies  by 
the  GPS  satellites.  The  receiver  produced  telemetry  inputs  to  the  DGD  measurement  unit  over  a 
voltage  range  (-5  to  +.')  volt.s)  th.it  was  twice  that  expected  (0  to  5  volts)  by  the  unit.  A  hardware 
modification  in  the  field  corrected  the  problem,  and  a  software  modification  pennitted 
recalibration  of  the  erroneous  data.  The  procedure  was  to  convert  the  erroneous  DGD  values 
back  to  voltages  and  then  to  recompute  the  correct  DGD  and  the  corresponding  TEC.  The  DGD 
unit  was  able  to  tolerate  some  negative  voltage,  although  not  as  much  as  -5  volts.  Thus,  the 
original  DGD  values  were  correct  except  at  certain  voltage  levels.  The  recalibrated  TEC  values 
still  were  somewhat  unreliable  near  those  levels,  but  we  were  able  to  refine  them  by  means  of  the 
DCP  measurements.  To  do  so,  we  obtained  reference  points  for  absolute  TEC  from  the  reliable 
portions  of  the  DGD  records  and  connected  them  by  means  of  the  (relative,  but  highly  precise) 
DCP  records. 
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III.  DATA  PROCESSING 


The  data  illustrated  in  Fig.  1  were  recorded  on  magnetic  tape  for  later  calibration 
refinement,  including  linearization,  and  further  processing.  Regarding  signal  fluctuations,  for 
instance,  the  intensity  .scintillation  inde.x  was  computed  as  follows: 

$4  =  <(P  -  <P>)2>>/2/<p2>l/2  ,  (1 ) 

where  <>  rcprc.scnts  averaging  over  1  min.  These  indices  are  .shown  in  the  top  and  bottom 
panels  of  Fig.  3  for,  respectively,  the  L2  and  LI  data  illustrated  in  Fig.  1.  The  faithful  manner  in 
which  the  two  indices  track,  that  for  the  lower  frequency  being  consistently  larger  than  that  for 
the  higher,  is  an  indication  of  high  data  quality.  The  obvious  departure  in  this  tracking  fidelity  at 
the  end  of  the  record,  when  the  receiver  was  about  to  be  switched  to  another  satellite,  gives  ready- 
indication  that  the  extreme  S4  value  for  L2  (well  in  excess  of  unity)  at  that  time  is  faulty.  The 
(valid)  L2  value  of  S4  approaching  0.9,  near  the  middle  of  the  time  span,  is  rare  for  L-band 
scintillation  at  high  latitudes,  although  it  is  encountered  more  frequently  near  the  geomagnetic 
equator,  especially  near  solar  maximum. 

While  scintillation  is  measured  at  each  of  the  two  frequencies.  TEC  is  determined  from 
the  DCP  and  DGD  between  them,  as  follovv.s.  Tlie  TEC  is  related  to  the  phase-path  defect.  0,,  on 
each  of  the  two  frequencies,  fj,  by 

TEC  =  Were,  (2) 

where  c  is  the  speed  of  light,  and  is  the  classical  electron  radius.  To  measure  the  signal  phase, 
one  needs  a  phase  reference  from  the  satellite.  To  provide  such  a  reference  is  the  reason  for 
receiving  both  LI  and  L2,  which  are  coherent  harmonics  of  a  common  oscillator.  Tlic  receiver 
produces  a  baseband  signal  who.se  phase  is  a  scaled  difference  between  <>1  and  Or  Tlie  two 
frequencies  are  sufficiently  clo.se  together  that  their  pha.se  variations  are  correlated  under  most 
mid-latitude  and  high  latitude  conditions.  In  this  circumstance,  the  pha.se  (in  radians;  obtained 
from  the  baseband  signal  and  output  from  the  receiver  is 

DCP  =  [l-(f2/f,)2j02.  (3) 

Combining  Equations  (2)  and  (3),  with  the  former  evaluated  at  f2.  yields 

TEC  =  0.37  DCP,  (4) 

in  TEC  units. 

The  DCP  is  free  of  effects  dtie  to  (the  time  integral  of)  geometrical  dopplcr  and  phase 
instability  in  the  transmitter  and  receiver  (other  than  differential  ones  between  the  two  carrier 
channels).  Its  time  history  yields,  via  Eq.  (4),  a  very  accurate  mca.surc  of  the  time  history  of 
TEC  variations  along  the  radio  propagation  path.  For  typical  values  of  TEC  (which  range 
globally  from  a  few  units  to  the  order  of  100  units),  however,  the  DCP  usually  is  ambiguous  by 
several  factors  of  2it.  TIic  ambiguity  is  resolved  by  means  of  the  DGD,  which,  when  measured 
in  nsec,  is  related  to  TEC  by 
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Figure  3.  Intensity  scintillation  indices,  S4,  computed  for  L2  (top)  and  LI  (bottom)  from  the  raw 
data  illustrated,  in  part,  in  the  bottom  two  panels  of  Fig.  1.  Time  span  is  13:13  through 
16;  13  UT.  The  three  spikes  in  S4  correspond  to  the  three  patches  of  scintillation  in  Fig. 
1.  The  large  spike  in  the  L2  value  occurred  later  than  the  end  of  Fig.  1  and  is  an  obvious 
artifact. 
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TEC  =  2.86  DGD 


(5) 


units.  Thus,  DGD  provides  only  a  very  coarse  measure  of  TEC,  but  it  is  not  ambiguous  unless 
the  TEC  level  reaches  279  units.  The  combination  of  DCP  and  DGD  yields  an  estimate  of 
absolute  TEC  and  a  precise  record  of  its  variations. 


The  TEC  values  recorded  by  the  foregoing  means  during  the  satellite  pass  that  yielded 
the  intensity  scintillation  Indices  shown  in  Fig.  3  are  presented  in  Fig.  4,  along  with  the  azimuth 
and  elevation  angles  of  the  satellite  as  viewed  from  Lerwick.  The  double  trace  at  the  top  shows 
the  TEC  measured  on  the  slant  path,  the  smoother  one  being  obtained  from  the  DGD  and  used  to 
set  the  absolute  level  of  the  more  finely  stmctured  one,  which  is  derived  from  the  DCP.  The 
latter  has  been  divided  by  the  secant  of  the  ray-path  incidence  angle  on  the  ionosphere  to  convert 
the  slant-path  TEC  to  an  equivalent  vertical  TEC,  which  is  shown  by  the  single  curve  plotted 
beneath  the  double  trace.  The  azimuth  and  elevation  are  shown  respectively  by  the  x  and  + 
traces. 


Figure  4.  TEC  values  corresponding  to  S4  records  shown  in  Fig.  3,  derived  from  the  DGD 
(smoother  of  the  two  overlapping  traces  at  the  top)  and  DCP  (more  finely  structured 
trace)  measurements  shown,  in  part,  in  Fig.  1.  Lower  solid  trace  is  equivalent  vertical 
TEC  derived  from  DCP  by  removing  effect  of  changing  slant  path-length.  X  and  + 
indicate  satellite  azimuth  and  elevation,  respectively.  Time  span  is  13:14  through  16:14 
UT. 


9 


Figure  5  contains  a  synopsis  of  TEC  data  obtained  at  Lerwick  over  the  24-hour  period  in 
which  the  GPS  track  that  produced  the  data  shown  in  Figs.  1,  3,  and  4  occurred,  along  with 
satellite-location  information.  The  solid  traces  show  the  equivalent  vertical  TEC  derived  from 
DCP  and  calibrated  by  means  of  DGD  during  eight  separate  satellite  tracks  of  the  type  illustrated 
in  Fig.  4.  The  sub-ionospheric  latitudes  and  longitudes  of  the  satellite  tracks  (location  of  point 
where  raypath  penetrates  350  km  altitude)  are  shown  respectively  by  the  x  and  +  traces. 
Individual  track  durations  range  from  about  an  hour  to  about  three  hours,  the  two-hour 
(approximately  13  to  15  hours  UT  in  Fig.  5)  track  shown  in  Figs.  3  and  4  being  of  typical 
duration. 


TIME  (U.T. ) 


Figure  5.  Synoptic  chart  of  equivalent  vertical  TEC  values  (solid  curves)  derived  from  DCP  and 
calibrated  by  means  of  DGD  from  measurements  at  Lerwick  over  the  24-hour  period  (31 
December  1989,  UT)  in  which  the  satellite  track  producing  the  data  in  Figs.  1,  3,  and  4 
occurred,  x  and  +  indicate  satellite  latitude  and  longitude,  respectively. 

IV.  RESULTS 

A  single  day's  data  from  a  single  station  contains  a  mix  of  infonnation  about  TEC  and 
scintillation  in  relation  to  several  variables  of  interest,  a  point  illustrated  for  TEC  in  Fig.  5  and  in 
Fig.  6.  In  both,  we  see  generally  higher  values  of  TEC  in  the  middle  UT  hours  (which  also  is 
midday  in  local  time  at  Lerwick)  than  at  night.  This  general  behavior,  more  consistently 
portrayed  in  Fig.  6  than  in  Fig.  5,  stems  from  the  dominance  of  solar  ultraviolet  radiation  in 
producing  the  ionosphere.  This  diurnal  variation  is  convolved  with  any  true  longitudinal 
variation  during  :i  single  track  as  the  satellite  drifts  eastward  (e.g.,  the  track  through  03  UT  in 
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Fig.  6)  or  westward  (track  through  18  UT).  As  the  satellite  drifts,  so  does  the  radio  ray-path  and 
its  intersection  with  the  ionosphere,  which  thus  drifts  to  later  or  earlier  times  of  day,  respectively. 
Similarly,  the  satellite's  northward  or  southward  motion  causes  the  intersection  to  move  in 
geographic  and  geomagnetic  latitude. 


Figure  6.  Synoptic  chart  of  equivalent  vertical  TEC  and  satellite  latitude  and  longitude  for  24- 
hour  period  of  20  October  1989,  UT. 


Superposed  on  the  drift  of  the  penetration  point  in  latitude  and  longitude  is  the  drift  of  the 
ionospheric  plasma  itself.  For  instance,  in  the  second  half  of  the  day  shown  in  Fig.  6, 
considerable  structure  may  be  seen  in  the  TEC  records.  This  very  likely  is  the  signature  of 
plasma-density  structures  convecting  in  the  dusk  sector  of  the  auroral  and  sub-auroral  F  layer. 
The  day  shown  (20  October  1989)  was  the  first  day  of  a  major  magnetic  storm  that  lasted 
through  the  next  day  and  on  into  23  October.  The  storm's  sudden  commencement  occurred  at 
09:17  UT  on  the  20th,  and  Kp  reached  8+  for  the  period  18  to  21  UT.  The  day  shown  in  Fig.  5 
(31  December  1989)  also  was  magnetically  disturbed,  especially  in  its  early  hours. 

Clearly,  detailed  analysis  of  TEC  and  scintillation  data  within  the  context  of  ionospheric 
dynamics  requires  separation  of  the  temporal  and  spatial  behaviors  described  above.  Such 
detailed  analysis  was  beyond  the  scope  of  this  preliminary  investigation,  but  we  were  able  to 
identify  events  of  interest  for  such  analysis  in  the  future.  For  example.  Fig.  7  contains  the  quick- 
look  spreadsheet  output  from  Lerwick  for  the  month  of  September  1989  (the  month  in  which 
observations  started).  We  note  four  days  in  which  substantial  signal  fluctuations,  characterized 
by  TEC  activity  at  level  3  or  above  and/or  by  scintillation,  took  place.  Three  of  these  days  --  the 
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Figure  7.  Quick-look  spreadsheet  output  from  Lerwick  for  September  1989,  indicating 
scintillation  and/or  elevated  TEC  activity  on  UT  dates  16,  18, 19,  22,  and  26,  all  of  which 
were  geomagnetically  disturbed  days. 
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18th,  19lh,  and  26th  -  were  the  three  most  disturbed  days  of  the  month  magnetically  (Kp  sum  = 
36+,  38-,  and  35-t-,  respectively).  The  remaining  one  -  the  22nd  -  was  the  next  most  disturbed 
day  magnetically  (Kp  sum  =  28+)  in  the  half  of  the  month  for  which  we  have  full  GPS  data  from 
Lerwick,  and  the  sixth  most  disturbed  in  the  month.  (The  fourth  most  disturbed  day  of  the 
month,  with  Kp  sum  -  3.5,  was  the  15lh.  The  greatest  magnetic  activity  on  that  day  occurred  in 
its  closing  hours,  when  the  GPS  receiver  was  not  operating,  as  indicated  by  the  Os  in  Fig.  7.  The 
fifth  most  disturbed  day  was  the  4th,  with  Kp  sum  =  29-,  prior  to  commencement  of 
observations.)  We  present  here  a  preliminary  discussion  of  signal  perturbations  on  22 
September,  during  which  scintillation  occurred  with  only  modest  TEC  activity  (level  1),  and  on 
26  September,  during  which  scintillation  occurred  with  TEC  activity  as  high  as  level  4. 

Figure  8  contains  the  synoptic  TEC  chart  from  22  September.  Although  classified  as 
having  only  level- 1  TEC  activity,  the  first  GPS  track  after  noon  UT  clearly  was  the  most 
disturbed  of  the  24-hour  period  shown.  It  was  during  the  latter  half  of  this  track  that  the 
scintillation  indicated  in  Fig.  7  occurred.  We  have  performed  further  analysis  of  the  data  from 
this  period.  The  latitude  and  longitude  tracks  (x  and  +,  respectively)  indicate  that  the  satellite 
was  moving  slowly  equatorward  and  eastward  across  the  arctic  circle  east  of  Iceland,  being  the 
highest-latitude  track  of  the  day.  The  ionospheric  penetration  point  scutined  over  the  Norwegian 
Sea,  in  the  general  vicinity  of  the  Faroe  Islands  northwest  of  the  Shetlands.  The  southernmost 
segment  of  this  latitude  span  had  been  sampled  somewhat  farther  to  the  east  (and,  therefore,  at  a 
slightly  lower  geomagnetic  latitude)  on  the  last  track  before  nooti  UT  (several  hours  earlier), 
with  generally  quiet  ionospheric  conditions  being  encountered  there  and  to  the  .south. 


Figure  8.  Synoptic  chart  of  equivalent  vertical  TEC  and  satellite  latitude  and  longitude  for  24- 
hour  period  of  22  September  1989,  UT.  Scintillation  occurred  during  UT  hours  14  and 
1 5,  in  the  pre.sence  of  level- 1  TEC  activity. 
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It  happens  that  one  of  the  "Transit"  satellites  of  the  Navy  Navigation  Satellite  System 
(NNSS)  was  tracked  from  Lerwick  near  the  time  of  both  of  the  foregoing  final  forenoon  and  fir^t 
afternoon  GPS  tracks  (L.  Kersley  and  G.  Bishop,  private  communication,  1991).  The  Transit 
receiver  output  the  dispersive  phase  (DCP)  between  15()-MHz  and  400-MHz  carriers,  which  are 
broadcast  from  all  the  NNSS  satellites.  Since  the  Transit  satellites  are  in  polar  (9()-deg 
inclination)  orbits  above  (1000  km)  the  ionospheric  F  layer,  their  dispersive-phase  records 
provide  a  latitudinal  scan  of  relative  TEC  in  the  vicinity  of  the  receiving  station.  By  means  noi 
known  to  us,  the  Transit  records  are  calibrated  to  an  aksoliite  level  of  TEC,  oin  interest  here 
being  only  in  TEC  gradients. 

The  latitudinal  TEC  scans  from  Transit  most  nearly  coincident  with  the  final  forenoon 
and  first  afternoon  GPS  tracks  on  22  September  1989  are  reproduced  in  Fig.  9.  We  note  that  the 
transition  region  between  the  midlatitude  (higher  TEC)  and  subauroral  (lower  TEC)  ionosphere 
is  characterized  by  a  steeper  gradient  in  the  afternoon  pass  than  in  the  morning  pass.  All  other 
conditions  being  equal,  the  afternoon  period  (steeper  gradient)  would  have  been  more 
susceptible  to  development  of  plasma  structures  at  transition  latitudes  by  means  of  convective- 
interchange  instabilities  than  would  the  morning  period. 


40  45  50  55  60  65  70  75  80 

LATITUDE  (deg) 


Figure 


9  Latitudinal  .sctins  of  'I'EC  ob.serx'cd  at  Lerwick  by  means  of  |)olar-orbiting  "Transit" 
satellites  between  0906  and  0918  UT  (shallow  curve;  ;tnd  between  1518  and  1529  UT 
(steeper  curve)  on  22  September  1939.  The  former  was  appro.\im:itely  along  the  353-deg 
meridian,  and  the  latter  approximately  along  the  4-deg  meridian.  (Data  provided  by  L. 
Kersley  via  G.  Bishop.) 
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The  afternoon  GPS  track  encompassed  a  latitude  span  some  ten  or  so  degrees  higher, 
geomagnetically,  than  did  tlie  morning  track.  Thus,  it  would  be  expected  generally  to  have 
sampled  a  more  dynamic  ionospheric  region.  On  22  September,  the  three-hour  planetary 
geomagnetic  activity  index,  Kp,  developed  as  foIlow.s:  '  ^  5-,  4+,  4+,  4-,  3-.  Thus,  the  last 

morning  GPS  track  occurred  in  a  developing  stage  (  •  ihrce-hour  period)  of  geomagnetic 

activity,  and  the  first  afternoon  pass  occurred  in  a  «!  voJoped  stage  (5ih  period).  By  the 
latter  time,  the  auronti  and  subauroral  regions  p.ohj-*..ib>  •.  ,<*  ,xpanded  eqiiatorward  *rom  their 
quie.scent,  quiet-time  locations. 

In  contrast  to  the  single  disturbed  GPS  track  on  j .  September,  all  post-noon  passes  on  26 
September  were  disturbed,  primarily  in  terms  of  TEC  fi..-  <  <ation.  Scintillation  was  encountered 
also.  The  persistence  of  the  post-noon  TEC  activity  i.  ^.iOent  in  Fig.  iO,  which  contains  the 
.synoptic  TEC  chart  for  the  day.  Noting  that  the  track  pat..’.rn  in  latituue  and  longitude  is  virtually 
identical  to  that  of  four  days  earlier  (Fig.  8),  we  see  that  the  TEC  structures  extended  more 
eqiiatorward  on  the  26th  than  on  the  22nd. 


Figure  10.  Synoptic  chart  of  equivalent  vertical  TEC  and  satellite  latitude  aiid  longitude  for  24- 
hour  period  of  26  September  1989,  UT. 

Thorough  investigation  of  signal  behavior  during  disturbed  conditions  would  employ  a 
combination  of  the  routinely  processed  data  and  the  raw  data,  recorded  at  20  sps.  Spectral 
analysis  of  the  raw  DCP  data  also  should  be  interpreta!’ 1?  in  terms  of  the  spatial  .spectrum  of 
relative  fine  structure  in  TEC  and  o*'  the  plasma-density  irregularities  responsible  for  it.  Such 
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spectral  analysis  has  been  performed  on  data  from  26  September  to  demonstrate  utility  of  the 
database  collected  under  this  ‘ontract  for  the  purpose.  Full  accounting  for  the  effects  of 
diffraction  would  be  somewhat  complicated  in  the  event  of  strong  scintillation  (S4  >  0.5,  say). 
Such  events  arc  rare  in  the  •jatabase  collected  under  this  contract,  .^o  interpretation  in  terms  of  the 
spatial  spectra  of  plasma  stru.,tures  should  be  feasible  in  future  research.  In  tn  s  prehminary 
;  ’vestigation  we  confined  our  efforts  to  the  routinely  processed  data. 

Figures  11  and  12  respectively  contain  the  TEC  and  scintillation-index  plots  for  the  first 
afternoon  pass  on  22  September.  VVe  see  that  scintillation  occurred  in  three  patches.  T'.;  first 
was  coincident  with  some  evident  fine  structure  n  TEC  and  the  second  with  a  local  TEC 
enhancement  en.  antered  near  1413  UT.  The  third  occurred  during  increasing  TEC  near  the  end 
of  the  pass.  The  presence  of  intensity  sc..iiiilation  means  that  the  individual  pl.ases  of  the  two 
GPS  signals  observed  at  the  ground  would  not  perfectly  track  the  variations  in  EC  imposed  by 
the  line-of-sight  integral  of  plasma  densit'  directly.  The  scintillation  is  sufficiently  weak, 
however,  that  the  effects  of  propagation  on  the  phases  of  signals  so  clo.sely  spaced  in  frequency 
probably  were  correlated,  a  conjecture  that  can  be  tested  by  means  of  the  individual  (raw) 
intensity  jecords.  If  so,  the  DCP  record  still  should  be  interpretable  in  terms  of  TEC,  and  its 
spectrum  in  terms  of  the  spatial  spectrum  of  plasma  structures. 


REF  VALUE: 49. 02 
PASS  FILE: 10 


SV  BIAS  CORP^CTION:  0.83 


Figi’re.  11,  TEC  measured  during  the  first  afternoon  GPS  track  in  Fig.  8. 
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SCFNTILLflTON  L2  SCINTiLLflTON  LI 


Figure  12.  Intensity  scintillation  indices,  S4,  measured  during  the  first  afternoon  GPS  track  in 
Fig.  8. 
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Equally  interesting  features  appear  in  the  I'EC  records  from  the  afternoon  of  26 
September,  when  the  line  of  sight  to  four  sequentially  observed  GPS  satellites  scanned  first  the 
highest-latitude  region  observable  from  Lerwick  and  then,  thrice,  lower  transitional  latitudes. 
The  TEC  record  from  the  first  (and  longest)  of  the  transitional  scans  is  shown  in  Fig.  13.  Of 
particular  interest  is  the  ten-minute-or-so  period  centered  on  about  17;27  UT,  when  strong  quasi¬ 
period  TEC  variations  were  recorded.  These  variations  could  be  the  signature  of  a  passing 
gravity  wave  launched  from  the  auroral  oval  under  geomagneiically  disturbed  conditions 
characterized  by  Kp  =  6+. 

REF  VRLUE:32.00 
PftSS  FILE:9 


SV  BIAS  CORRECTION:  0.00 
o 


Figure  13  TEC  measured  during  the  second  afternoon  GPS  track  in  Fiii.  10.  Note  strong  quasi- 
periodic  signature  centered  on  about  17:27  UT. 


V.  CONCLUSION  AND  DISCUSSION 

In  this  preliminary  investigation,  we  have  found  decidedly  more  variation  in  DCP,  DGD, 
and  L-band  signal  intensity  in  the  database  from  the  Shetland  Islands  than  in  that  from  Hanscom 
Field  This  behavior  di.scloses  a  far  greater  incidence  of  sharp  TEC  gradients  and  scintillation 
resulting  therefrom  in  the  ionospheric  region  viewed  from  Lerwick  than  in  that  viewed  from 
Hanscom.  New  studies  that  have  been  propo.sed  -  one  employing  routinely  proces.sed  TEC  and 
scintillation  data  of  the  type  displayed  in  this  report  and  one  employing  the  raw  (20-sp.s)  data 
collected  under  this  contract  -  should  disclose  the  nature  of  the  ionospheric  proccs.ses 
responsible  for  this  observed  difference. 


The  fields  of  view  of  the  aforesaid  two  databtises  are  well  separated  in  longitude  and 
largely  overlapping  in  latitude.  Still,  their  center  points  are  separated  in  geomagnetic  latitude  by 
about  five  degrees,  and  this  could  be  enough  to  account  for  the  substantially  greater  degree  of 
activity  seen  at  Lerwick  than  at  Hanscom.  To  test  this  possibility  and  to  assess  any  residual 
effect  due  to  longitude,  each  routinely  processed  data  segment  needs  to  be  placed  in  a 
location/time  bin.  The  bin  should  be  defined  by  the  latitude  and  longitude  of  the  penetration 
point  of  the  radio  line  of  sight  in  the  F  layer  (at  350  km,  say)  and  the  local  time  at  that  point. 
This  is  a  simple  matter  of  geometiy.  The  local  time  and  season  then  should  be  used  to  locate  the 
penetration  point  relative  to  the  sohir  terminator,  so  that  gradients  associated  therewith  can  be 
identified. 

Each  location/time  would  correspond  to  a  particular  geomagnetic  latitude  and 
longitude/time,  which  also  should  be  identified  using  a  standard  model  of  the  geomagnetic  field 
(preferably,  the  International  Geomagnetic  Reference  Field).  Given  the  clear  evidence  in  the 
Lerwick  database  of  geomagnetic  control  over  the  TEC  and  scintillation  activity  observed  there, 
the  data  segments  also  should  be  placed  into  Kp  bins.  It  may  be  that  another  inde.x  of 
geomagnetic  activity,  such  as  the  auroral  electrojet  index,  woidd  parameterize  the  activity  more 
systematically.  This  preliminary  analysis  disclosed  an  unusually  clear  relationship  between  the 
activity  and  Kp,  however,  so  the  latter  would  provide  a  useful  starting  point. 

Detailed  analysis  of  selected  portions  of  the  raw  data  would  provide  additional 
infonnation  on  signal  behavior  and  the  underlying  plasma  structures  that  cau.se  it.  As  a  first  step, 
the  selected  intensity  and  DCP  measurements  should  be  plotted  for  visual  inspection  at  full  (20- 
sps)  resolution.  Such  inspection  of  the  DCP  records  may  yield  identifiable  signatures  of  TEC 
fine  structure  due  to  particular  dynamical  processes,  analogous  to  that  identified,  at  coarser 
resolution,  in  Fig.  13  as  possibly  due  to  a  passing  gravity  wave.  Inspection  of  the  corresponding 
intensity  records  would  disclo.se  propagation  effects  that  might  have  to  be  accounted  for  in 
interpretations.  Useful  in  this  regard  would  be  in.spection  of  both  the  LI  and  L2  records, 
performing  a  "correlation  by  inspection"  to  ferret  out  any  instances  in  which  DCP  may  be  altered 
by  decorrelation  of  propagation  effects  on  the  two  frequencies.  We  believe  that  such  instances 
have  been  rare.  Whether  such  propagation  effects  are  ab.sent  or  are  to  be  dealt  with,  the  DCP 
records  should  be  interpreted  either  deterministically  or  statistically,  depending  upon  the 
character  of  a  particular  event.  In  the  latter  case,  spectral  analysis  should  be  u.sed  to  interpret  the 
DCP  records  in  terms  of  the  spatial  spectrum  of  plasma-density  structures,  and  we  have  proven 
its  feasibility  with  the  database  collected  under  this  contract. 

It  is  recommended  that  the  databa.se  dc.scribcd  herein  be  incorporated  into  a  more 
extensive  one  containing  TEC  and  .scintillation  data  from  Han.seom  Field;  the  Shetland  Islands; 
and  Thule,  Greenland.  We  further  recommend  analysis  t)f  the  contents  of  the  combined  database 
to  fomi  conclusions  about  the  solar-maximum  behavior  of  TEC  and  its  fluctuations  at  trough- 
edge,  auroral,  and  polar-cap  latitudes.  The  results  of  such  an  analysis  would  contribute  to 
software  intended  to  characterize  the  transionospheric  radio  propagation  channel  employed  for 
communication  and  sur\'cillancc  purpo.scs. 
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